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ABSTRACT 

We present high precision radial velocity observations of HD 17156 during a transit of its eccen- 
tric Jovian planet. In these data, we detect the Rossiter-McLaughlin effect, which is an apparent 
perturbation in the velocity of the star due to the progressive occultation of part of the rotating stel- 
lar ph otosphere by the transiting planet. This system had previously been reported by iNarita et al.l 
(2008) to exhibit a A = 62° ± 25° misalignment of the projected planetary orbital axis and the stellar 
rotation axis. We model our data, along with the INarita et alJ data, and obtain A = 9.4° ± 9.3° for 
the combined data set. We thus conclude that the planetary orbital axis is actually very well aligned 
with the stellar rotation axis. 

Subject headings: stars: individual(HD 17156) — planetary systems 



1. INTRODUCTION 

Transiting extrasolar planets allow us to perform crit- 
ical tests of models of planetary system formation and 
evolution. In our solar system, all of the planets or- 
bit approxima tely, but not exactl y, in the solar equa- 
torial plane. iBeck fc Giles! (|2005f) find that the angle 
between the plane of the ecliptic and the solar equator 
is 7.155 ± 0.002°. The near co-planarity of the plan- 
etary orbits in our solar system has influe nced models 
of p lanetary system origin from the time of iKantl (|1755f ) 
and Laplace (1796) to essentially all modern models (e.g. 
lLissauerlll995HPollack et al .11 1 9961 : [Bossll2000h . Given the 
level of misalignment in our solar system, and the obser- 
vation of war ps in debris disks around nearby stars such 
as Pic (e.g. iBurrows et al.l ll995t iMouillet et allll997t 
iHeap et all 12000). it is obvious that there are common 
processes which give rise to some small level of spin-orbit 
misalignment in planetary systems. The degree of mis- 
alignment in real planetary systems must depend on the 
initial conditions (the initial asymmetries of the collaps- 
ing cloud), the physics of disk formation and evolution, 
and the physics of stellar mass and angular momentum 
loss during the T Tauri phase. Most transiting planets 
have periods of just a few days and orbits of low eccen- 
tricit y. Significan t exce ptions are HD 147506b (HAT-P- 
2b) (|Bakos et al.l 120071 ) . with a 5.6 day orbital p eriod 
and eccentricity of 0.52, HD 17156 ([Fischer et all 120071 ; 
iBarbieri et al.l l2007h . which has the longest orbital pe- 
riod (21 days) and the largest eccentricity (0.67) of all of 
the transiting exoplanets, and XO-3b, a massive (13.25 
Mj) planet which has an eccentricity of 0.26 in spite of 
it short orbital period of 3.192 days ( Johns-Krul l et al.l 
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2008) . Scenarios for the formation of short-period plan- 
etary systems do not necessarily deliver those planets 
to their present semi-major axes with low eccentricity. 
Type II migra tion can result in plan ets with moderate 
eccentricities ([Sari fc Goldreichll2004h . while dynamical 
intera ctions among newly form ed planets and planetes- 
imals (|juric &; Tremainej|2007f) or the Ko zai mechanism 
([Holman et all 119971 : IWu fc Murray! 12003ft can result in 
large eccentricities and significant inclination of plane- 
tary orbital planes from the plane of the stellar equa- 
tor (which is presumably close to the plane of the inner 
portion of the proto-planetary disk). While the shortest 
period transiting systems have probably been tidally cir- 
cularized, longer period systems can easily retain large 
eccentricities over the main-sequence lifetime of the par- 
ent star. A planet that has undergone significant gravita- 
tional scattering or Kozai excitation would not necessar- 
ily retain a low inclination relative to the stellar equator. 
Thus, the report of a possible spin-o rbit misalignment o f 
the HD 17156b transiting planet bv INarita etaLI (|2008l ) 
is extremely interesting, and calls for further detailed 
investigation. In this paper, we report our own spectro- 
scopic observations of a transit of HD 17156 by its planet, 
using two telescopes at McDonald Observatory. 

2. OBSERVATIONS 

We observed the transit of HD 17156 by its hot-Jupiter 
companion on the night of 25 December 2007 UT, us- 
ing both the 2.7m Harlan J. Smith Telescope (HJST) 
and the Hobby-Eberly Telescope (HET) at McDonald 
Observatory. The HJST obser vations used the 2dcoude 
spectrograph ([Tull et alj|1995| ) in its "F3" mode, which 
gives a spectral resolving power of R = X/SX = 60,000. 
This mode is referred to as "cs23". A temperature sta- 
bilized I2 gas absorption cell is used to impose the ve- 
locity metric for precise radial velocity measurements of 
the stellar spectrum. Details of the 2.7m cs2 3 observing 
and data re duction procedures are given bv lEndl et "all 
(|2004l [20061) . Velocity observations were started before 
the expected beginning of the transit, and were contin- 
ued until after the expecte d end of the transit. A t total 
e $M#tef^te^a#fl x !S^nmte^^ 

Table [T] gives the relative radial velocities for the HJST 
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TABLE 1 
McDonald Observatory 2.7m HJST 
Relative Velocities for HD 17156 
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TABLE 2 
Hobby-Eberly Telescope HRS 
Relative Velocities for HD 17156 
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observations of HD 17156. 

HET observations were made on the same night 
(25 Dece mber 2007 U T) using the High Resolution Spec- 
trograph (|Tull 119981 ) in its R = 60,000 mode. Due 
to its fixed-zenith-distance design, we were only able to 
observe HD 17156 from shortly before the beginning of 
the transit to just past past mid-transit. The observa- 
tions were planned to obtain as many 600 s exposures of 
HD 17156 as possible during the 2.1 hour track length. 
The 13 th target exposure on the HET was terminated af- 
ter 455 s when the fiber-instrument-feed reached the end 
of track. Details of the instrument configuration and 
the data reduc t ion and ana lysis procedures are given by 
ICochran et ail ((20041 [2007T) . Table CU gives the relative 
radial velocities for the HET observations of HD 17156. 

For both the HJST and HET data, observation times 
and velocities have been corrected to the solar system 
barycenter. The uncertainty a for each velocity in the 
table is an internal error computed from the variance 
about the mean of the velocities from each of the ~ 2A 
small chunks into which the spectrum is divided for the 
velocity computation. Thus, it represents the relative 



uncertainty of one velocity measurement with respect to 
the others for that instrument, based on the quality and 
observing conditions of the spectrum. This uncertainty 
does not include other intrinsic stellar sources of uncer- 
tainty, nor any unidentified sources of systematic errors. 
The two different spectrographs have independent arbi- 
trary velocity zero points, and thus there is some con- 
stant offset velocity (determined below and denoted as 
7) between the data sets presented in Tables Q] and [2] 

3. ROSSITER-MCLAUGHLIN EFFECT MODEL 

A variety of different types of models have been 
used by others to ana lyze observat i ons of the Rossiter- 
McLaughlin (RM) (jRossiterl 119241: iMcLaughlinl Il924f ) 
effect for transiting planets. iQueloz et al.l (|2000f ) di- 
vided a model stellar photosphere into a large num- 
ber of cells, and then used a "Gaussian shape cross- 
correlation model" with a linear lim b darkening law to 
compute the radial velocity anomaly. lOhta et "aLl (|2005f ) 
developed analytic expressions for the apparent radial 
velocity perturbatio n during a trans it, in several differ- 
ent approximations. iGimenezl (|2006[ ) developed another 
set of analytic expressions for the RM effect which uti- 
lize a more generalized higher order limb darkening ex- 
pression. A more elaborate technique was developed by 
W inn et all (|2005f ). who first computed an approxima- 
tion to the disk-integrated stellar spectrum. They then 
computed a Doppler shifted and intensity scaled spec- 
trum of the portion of the disk that would be blocked 
by the transiting planet and subtracted this from their 
disk-integrated spectrum. This spectrum was then mul- 
tiplied by their high-resolution iodine spectrum, and the 
result was processed through their radial velocity code 
to compute model velocities in the same manner as the 
observed data. 

We analyzed our data using a model that is a hybrid of 
these methods. We start ed by adopti n g the HD 17156b 
system parameters from iNarita et al.1 (1200 8) . We then 
computed the orbit of the planet around the star, as 
we would view it from Earth. This gave us the appar- 
ent offset of the planet from the center of the star as 
a function of time through the transit. We divided the 
st ellar disk into a 400 x 400 grid of cel ls, in the manner 
ofjOrteloz etaLl (l200l . iSnellenl (l200l . or iWinn et all 
(2005). For each photospheric cell, we computed a spe- 
cific intensity usi ng the non-li near four-parameter limb 
darkening law of IClaretl (|2000T ). Each photospheric cell 
is also assigned a radial velocity due to both the stellar 
orbital motion and the stellar rotation, with the stel- 
lar -usini as a model parameter. For each time step 
during the transit, from first contact to fourth contact, 
we compute which stellar photospheric cells are blocked 
by the transiting planet. We then integrate the un- 
blocked Doppler-shifted and intensity weighted stellar 
photospheric cells to compute both the RM radial veloc- 
ity perturbation during the transit and the transit pho- 
tometric lightcurve. 

We fully recognize the limitations and approximations 
inherent in this modeling procedure. First, there is no 
a priori reason to assume that limb-darkening should 
follow any particular law. Also, the limb darkening in 
photospheric absorption lines is quite different from the 
limb darkening in the con tinuum . Whil e this appears to 
be taken into account by IClareH (|2000f ). the limb dark- 
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TABLE 3 

Assumed System Parameters 



Parameter 



Value 



Source 



M* (M ) 
R* (Re) 

Claret al 
a2 
a3 
a4 

Rp (Rjup) 

P (days) 
T p (HJD) 
K (ms" 1 ) 

e 

lu (degrees) 



1.2 
1.47 
0.5346 
0.1041 
0.4189 
-0.2584 
1.01 
21.21691 
2453738.605 
273.8 

0.670 
121.3 
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ening parameterization is based on model atmospheres 
rather than on real stars. More importantly, in our 
model we compute the specific-intensity weighted ap- 
parent Doppler shift of the visible portion of the pho- 
tosphere. On the other hand, our spectra record transit- 
perturbed stellar absorption line profile shapes from 
which we measure an apparent Doppler shift using a 
computer code that assumes an unperturbed line-profile 
shape. In future improvements to our RM model, we will 
attempt to improve several of these limitations. 

4. DATA ANALYSIS 

We used the model described in Sec. [3] to analyze si- 
multaneously the data sets from the HJST c s23, the HET 
HRS, and the observations published by iNarita et al.l 
( 2008), which we will refer to as the "OAO/HIDES" data 
set. Since each data set has its own independent velocity 
zero-point, we allowed the systemic velocity of each data 
set to be an independent free parameter in the analysis. 
The values of the fixed parameters for the analysis are 
given in Ta ble [3l The planet ary orbital elements were 
taken from llrwin et al.l (|2008f ) . These elements are es- 
sentially indistin guishable from the single-planet fit of 
IShort et alj (|2008f ). We note that the conclusions of this 
work depend on the particular values of these parameters 
we adopt. If any of these parameters turn out to be in 
significant error, those errors will propagate through this 
analysis. 

In modeling the data, we allowed the orbital plane in- 
clination z, the projected angle A between the planetary 
orbital axis and the stellar rotation axis, the projected 
stellar rotational velocity vsmi, as well as the systemic 
velocity of each separate data set to be free parameters. 
We then minimized the chi-squared of the model fit to 
the data. 

For the combined data sets, we obtained A = 9.4°±9.3° 
at vsmi = 6.3 ± l.lkms^ 1 and i = 85.9° ± 0.4°. This 
indicates that the planet is orbiting near the stellar equa- 
torial plane, to within the uncertainty of our determina- 
tion. Our derived usini i s somewhat larger than the 
value of 4.7 ms -1 fo und by INarita et al (|2008t ) and the 
2.6 km s given bv lFischer et all (l2007n . Our inclina- 
tion is within one a of the INarita et al.l value of 85.65°. 
However, our value for A from the combined data differs 
significantly from that of lNarita et alJ (|2008h , who found 
A = 62° ±25°. 

In order to understand the reason for this difference, 
we then modeled each of the three data sets separately. 




0.98 0.985 0.99 

Orbital Phase 

Fig. 1. — The fit of our Rossiter-McLaughlin effect model (solid 
black line) to the observational data sets. The HJST cs23 data are 
shown as the (red) Xs, the HET HRS data are (blue) solid squares, 
and the OAO/HIDES data are (green) open triangles. The dashed 
black line is the orbital velocity of the star in the absence of any 
Rossiter-McLaughlin velocity perturbation. 

The derived A and i for the combined fit, as well as for 
each individual data set are given in Tab. |U We also 
give 7c, the systemic velocity for each data set in the 
combined fit, as well as 7j, the systemic velocity for each 
individual data set fit. From the OAO/HIDES data set 
alone, we computed A = 59.4° ± 19.7° at i = 86.2° ± 
0. 8°. If we fix the in clination at the i = 85.65° value 
of lNarita et all (pool ), we get A = 61.4° ± 28.5°. The 
surprisingly excellent agreement of all of these value for 
the OAO/ HIDES data set thus validates the modeling 
process of Na rita et al.l . We note that the results from 
the HJST data alone are in very good agreement with 
the combined results. Due to the design limitations of 
the HET, the HET/HRS data only covered the first half 
of the transit. Thus, the code could trade-off A vs. the 
systemic velocity for the data set, and derived a slightly 
lower x 2 with A = -32.4°±25.2° for the HET data alone. 

The model fit to the data is shown in Fig. [T] We 
also computed a model with no RM effect by setting 
vsmi = 0. This removes the stellar rotation, and thus 
there is no apparent Doppler shift of the portion of the 
stellar photosphere eclipsed by the planet. The v sin i = 
model gave \ 2 = 71.56 (Xr = 1-43) as opposed to 
X 2 = 40.53 (Xr = 0-83) for the best fitting RM model for 
the combined data set. Thus, we conclude that the RM 
effect was indeed convincingly detected in spite of the ap- 
parent noise in the data. Examining the x 2 of each indi- 
vidual data set in the full model and the v sin i = model 
showed that even in the noisiest case of the OAO /HIDES 
data, the RM effect reduced the total % 2 by 13.1. 

Due to the symmetric signature of the RM effect in the 
A«0 case, it is critically important to sample the entire 
transit. This is evidenced by the fit to the HET/HRS 
data alone, in which despite the high quality of the data, 
only the first half of the transit was observable, and 
therefore erroneous parameters are derived. In order to 
successfully calibrate the baseline radial velocity varia- 
tion, observations well before and after the transit are 
required. If the radial velocity can be well calibrated 
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TABLE 4 

Rossiter-McLaughlin Model Results 



Data Set 


A 


i 


vs'mi 


x 2 


degrees of 


7c 


7< 




(degrees) 


(degrees) 


(kms -1 ) 




freedom 


m s — 1 


ms — 1 


Combined 


9.4±9.3 


85.9±0.4 


6.3±1.1 


40.53 


49 






HJST cs23 


4.5±15.6 


86.1±0.6 


7.1±2.1 


13.19 


11 


-7853.0 


-7854.0 


HET HRS 


-32.4±25.2 


86.3±0.5 


4.8±1.5 


3.71 


9 


73.3 


65.6 


OAO/HIDES 


59.4±19.7 


86.2±0.8 


8.6±1.9 


14.25 


21 


142.1 


146.5 



by observations before and after transit, the influence of 
uncertainty in the zero-point velocity offset (7) will be 
negligible on the RM effect parameters. 

5. DISCUSSION 

We have reanalyzed the data of iNarita et al.l ([2008) 
along with our own independent observations of the spec- 
troscopic transit of the eccentric exoplanet HD 17156b, 
and we find that A = 9.4° ± 9.3°. We conclude that this 
exoplanetary system is similar to almost all of the other 
short-period transiting exoplanetary systems studied us- 
ing the RM effect so far, in that it shows that the pro- 
jected planetary orbital axis appears to be aligned with 
the projected stellar rotation axis, to within our measure- 
ment precision. The only remai ning notable exceptio n is 
XO-3b, which was reported by iHebrard et all (|2008f ) to 
show a very significant misalignment of A = 70° ± 15°. 

The HD 17156b system is extremely interesting be- 
cause it has a very high eccentricity (e = 0.67) for such 
a short orbital period of 21.2 days. Most other transiting 
systems have much shorter orbital periods and eccentric- 
ities near zero. Many of the interesting ways of getting 
a planet into this type of orbit might well result in the 
possibility of an orbital plane significantly inclined to the 
stellar equatorial plane. Dynamical interactions with an- 
other nearby planet could result in ejection of the other 
body and a large orbital inclination of the surviving body. 
The pre sence of a th i rd pla net in the system, as sug- 



:ni rd p 

gested by IShort et ail (|2008h . would have some effect 



on the orb i tal el ements of HD 17156b. However, as 
I Short et"aT1 (|2008f ) discussed, the primary effects are a 
small oscillation of the eccentricity and a secular advance 
of the argument of periastron. As long as the two plan- 
ets are approximately coplanar, there would be no in- 
duced change in the orbital inclination. Thus, it appears 
that the dynamical process that placed HD 17156b into 
a short-period highly eccentric orbit did not significantly 
affect the orbital inclination of the system with respect 
to the stellar equatorial plane. 

We are extremely grateful for receiving Director's Dis- 
cretionary time on the HET on very short notice in 
order to obtain the data presented here. S.R. would 
like to acknowledge support provided by NASA through 
Hubble Fellowship grant HST-HF-01190.01 awarded by 
STScI, which is operated by AURA Inc., for NASA, un- 
der contract NAS 5-26555. This material is based on 
work supported by NASA under Grant NNG05G107G 
issued through the TPF Foundation Science Program, 
and under Cooperative Agreement NNA06CA98A issued 
through the Kepler Program. The Hobby-Eberly Tele- 
scope (HET) is a joint project of the University of Texas 
at Austin, the Pennsylvania State University, Stanford 
University, Ludwig Maximilians Universitat Miinchcn, 
and Georg August Universitat Gottingen. The HET is 
named in honor of its principal benefactors, William P. 
Hobby and Robert E. Ebcrly. 

Facilities: HET, McD:2.7m. 



REFERENCES 



Bakos, G. A., et al. 2007, ApJ, 670, 826 
Barbieri, M., et al. 2007, A&A, 476, L13 
Beck, J. C, & Giles, P. 2005, ApJ, 621, L153 
Boss, A. P. 2000, ApJ, 536, L101 

Burrows, J. C, Krist, J. E., Stapelfeldt, K. R., & WFPC2 

Investigation Definition Team. 1995, BAAS, 27, 1329 
Claret, A. 2000, A&A, 363, 1081 
Cochran, W. D., et al. 2004, ApJ, 611, L133 

Cochran, W. D., Endl, M., Wittenmyer, R. A., & Bean, J. L. 2007, 
ApJ, 665, 1407 

Endl, M., Cochran, W. D., Wittenmyer, R. A., & Hatzes, A. P. 

2006, AJ, 131, 3131 
Endl, M., Hatzes, A. P., Cochran, W. D., McArthur, B., Prieto, 

C. A., Paulson, D. B., Guenther, E., & Bedalov, A. 2004, ApJ, 

611, 1121 

Fischer, D. A., et al. 2007, ApJ, 669, 1336 
Gimenez, A. 2006, ApJ, 650, 408 

Heap, S. R., Lindler, D. J., Lanz, T. M., Cornett, R. H., Hubeny, 

I., Maran, S. P., & B, W. 2000, ApJ, 53 9, 435 
Hebrard, G., et al. 2008, arXiv:0806.0719]»>'l 

Holman, M., Touma, J., fe Tremaine, S. 1997, Nature, 386, 254 
Irwin, J., et al. 2008. larXiv: 0801.1496 
Johns-Krull, et al. 2008, ApJ, 677, 657 

Juric, M., & Tremaine, S. 2007, astro-ph/0703160, ApJ, in press 
Kant, I. 1755, General History of Nature and Theory of the Heavens 
(Koningsberg: Perersen) 



Laplace, P. S. 1796, Exposition du systeme du monde (Paris: 

Cercle-Social) 
Lissauer, J. J. 1995, Icarus, 114, 217 
McLaughlin, D. B. 1924, ApJ, 60, 22 

Mouillet, D., Larwood, J. D., Papaloizou, J. C. B., & Lagrange, 

A. M. 1997, MNRAS, 292, 896 
Narita, N., Sato, B., Ohshima, O., & Winn, J. N. 2008, Publ. 

Astron. Soc. Japan, 60, LI 
Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118 
Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 

Podolak, M., & Greenzweig, Y. 1996, Icarus, 124, 62 
Queloz, D., Eggenberger, A., Mayor, M., Perrier, C, Beuzit, J. L., 

Naef, D., Sivan, J. P., & Udry, S. 2000, A&A, 359, L13 
Rossiter, R. A. 1924, ApJ, 60, 15 
Sari, R., & Goldreich, P. 2004, ApJ, 606, L77 

S hort, P., W elsh, W. F., Orosz, J. A., & Windmiller, G. 2008, 

ErXiv:0803.2935] 
Snellen, I. A. G. 2004, MNRAS, 353, LI 

Toll, R. G. 1998, Proc. Soc. Photo-opt. Inst. Eng., 3355, 387 
Tull, R. G., MacQueen, P. J., Sneden, C, & Lambert, D. L. 1995, 

PASP, 107, 251 
Winn, J. N., et al. 2005, ApJ, 631, 1215 
Wu, Y., & Murray, N. 2003, ApJ, 589, 605 



